ASCA Observations of the Absorption-Line Features from the 
Super-Luminal Jet Source GRS 1915+105 
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ABSTRACT 



We have carried out a precise energy spectral analysis of the super-luminal 
jet source GRS1915+105 observed with ASCA six times from 1994 to 1999. The 
source was so bright that most SIS data suffered from event pileup. We have 
developed a new technique to circumvent the pileup effect, which enabled us to 
study the spectrum in detail and at high resolution (AE/E rs 2 %). In the 

m . 

Ka, Fe xxvi Ka, as well as blends of the absorption lines of Ni xxvn Ka + 
Fe xxv Kf3 and Ni xxvm Ka + Fe xxvi K/3, were observed. Such absorption 
lines have not been found in other objects, except for iron absorptions lines 
from GRO J1655— 40, another super-luminal jet source (Ueda et al. 1998). We 
carried out a "curve of growth" analysis for the absorption lines, and estimated 
column densities of the absorbing ions. We found that a plasma of moderate 
temperature (0.1-10 keV) and cosmic abundance cannot account for the observed 
large equivalent widths. The hydrogen column density of such plasma would 
be so high that the optical depth of Thomson scattering would be too thick 
(Ng > 10 24 cm~ 2 ). We require either a very high kinetic temperature of the ions 
(> 100 keV) or extreme over-abundances (> 100 Z Q ). In the former case, the 
ion column densities have reasonable values of 10 17 — 10 18 cm -2 . We modeled 
the absorber as a photo-ionized disk which envelops the central X-ray source. 
Using a photo-ionization calculation code, we constrain physical parameters of 
the plasma disk, such as the ionization parameter, radius, and density. Estimated 
parameters were found to be consistent with those of a radiation-driven disk wind. 
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These absorption-line features may be peculiar to super-luminal jet sources and 
related to the jet formation mechanism. Alternatively, they may be common 
characteristics of super-critical edge-on systems. 

Subject headings: binaries: close — stars: individual (GRS 1915+105) — X-rays: 
stars 



1. Introduction 

GRS 1915+105 was discovered as a transient X-ray source with WATCH/Granat (Castro- 
Tirado, Brandt, & Lund 1992), and later recognized as a super-luminal Galactic jet source at 
a distance of 12.5 kpc (Mirabel & Rodriguez 1994). The velocity of the jet and its inclination 
were estimated by Mirabel & Rodriguez (1996) to be 0.92 c and 70°, respectively . The opti- 
cal type of the companion star and binary parameters are as yet unknown. The formation of 
the relativistic jet is considered to be related to instabilities in a super-critical accretion disk 
(Belloni et al. 1997; Mirabel et al. 1998), but the precise formation mechanism is not yet 
clear. GRS 1915+105 is extremely important because it is a potential key to understanding 
the production of jets, which is relevant to astrophysical systems ranging from binaries to 
AGNs. To obtain information on the environment of the central engine of jet sources, X-ray 
spectroscopy with a fine resolution is desirable. 

GRS 1915+105 was observed with the X-ray astronomy satellite ASCA, and surprisingly, 
absorption features of calcium and iron K structure were found in the spectrum (Ebisawa 
1997a; Ebisawa 1997b; Kotani et al. 1997a; Ebisawa et al. 1998). These features were 
interpreted as absorption lines due to Ka-resonant scattering of helium-like and hydrogen- 
like ions (Kotani et al. 1999a; Kotani et al. 1999b). Such X-ray absorption lines have never 
been found in other X-ray binaries except for GRO J1655— 40, another super-luminal jet 
source (Ueda et al. 1997; Ueda et al. 1998). Ueda et al. (1997; 1998) identified helium-like and 
hydrogen-like iron-Ka absorption lines, and estimated the physical condition of the plasma 
scattering iron-Ka photons. Since two of the super-luminal jet sources show absorption lines, 
it is natural to consider the lines in the context of jet production. The absorption lines may 
be a byproduct of the unknown jet formation mechanism. In this paper, we discuss the 
identification and origin of the X-ray absorption lines in GRS 1915+105 in detail. 
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2. Observations 

GRS 1915+105 was observed with ASCA six times from 1994 to 1999 (Table 1). In this 
paper we mainly report and discuss the first two observations in 1994 and 1995, in which 
characteristic absorption lines were detected. In 1994, the source was observed as a Target- 
Of-Opportunity (TOO), responding to the report of a hard X-ray outburst (Sazonov et al. 
1994). The SIS was always operated in 1-CCD BRIGHT mode except for the observation in 
1995, when the SIS was in 1-CCD FAINT mode. The telemetry limit is 256 counts/s/SIS 
in BRIGHT mode at the highest telemetry rate, and 64 counts/s/SIS in FAINT mode. 
Therefore, the telemetry was saturated and a significant fraction of SIS events was lost in 
the 1995 observation. The GIS was always operated in PH mode with the nominal bit 
assignment (Ohashi et al. 1996). The source flux in the 2-10 keV band was 0.4, 0.8 and 0.9 
Crab for the observations in 1994 (Nagase et al. 1994), 1995 and 1996, respectively. Due 
to the high counting rate, the GIS spectra suffered from the both dead-time effects and 
telemetry saturation, which we correct for in the present analysis (Makishima et al. 1996). 

3. Data Reduction 

3.1. Extracting SIS and GIS Spectra 

Event pileup effects must be removed from SIS data to perform precise spectroscopy for 
bright sources like GRS1915+105. We have developed a new method which can correct both 
of the two symptoms of event pileup, namely, the reduction in counting rate and hardening of 
the spectrum. The correction method we applied is described in detail in Appendix A. The 
spectra taken in 1994 and 1996 were corrected with this code. As for the 1995 data, pileup 
correction is not required, because a significant fraction of the image including the part where 
pileup effects would be most severe was lost due to the telemetry saturation. The response 
function of the X-ray telescope (XRT) (Serlemitsos et al. 1996) for each selected region was 
calculated. Background subtraction was not performed since the source was brighter than 
the background by two orders of magnitude over the entire energy range. The gains of 
SIS-0 and SIS-1 were calibrated based on the instrumental gold-M-edge and silicon-K-edge 
features. The two spectra were then combined and a systematic error of 2 % was added to 
each spectral bin, in quadrature with the statistical error. 

GIS events were collected within a circular region with a radius of 6' centered on the 
source position for each sensor, and the corresponding XRT response functions were calcu- 
lated. Background was again not subtracted. After gain calibration using the instrumental 
gold M edge, the GIS-2 and GIS-3 spectra were combined. A systematic error of 2 % was 
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included for each energy bin of the spectra. 

The resultant SIS and GIS spectra were found to be not completely consistent with 
each other: The SIS spectra were systematically harder than the GIS spectra even after the 
pileup correction was applied. We also tried other pileup correction codes (e.g., Ebisawa 
et al. 1996; Ueda et al. 1997), but the results were essentially the same. Here we consider 
the GIS spectrum to be more reliable than the SIS to study the continuum spectral shape, 
since the SIS data for bright sources have calibration uncertainties other than pileup, e.g., 
uncorrectable DFE and the echo effect (Yamashita et al. 1997 and references therein). On 
the other hand, the gain dependence on count rate of the GIS is well understood (Makishima 
et al. 1996). We do not pursue the GIS and SIS spectral difference further in this paper, but 
rather discuss the spectral features common to both the SIS and GIS. Note that pileup is 
a global effect, so it will hardly affect any sharp spectral features like the absorption lines 
which we are mostly interested. 



4. Spectral Analysis Results 

4.1. Continuum Spectrum 

We tried several spectral models, and found that both the SIS and GIS spectra are 
roughly expressed with a partially attenuated cut-off-power-law model: 

exp[-<7(£) N m ] x {/ exp[-<7(£) N m ] + (1 - /)} x A exp[-E/E Md ] E~ T . 

Where E is the incident photon energy, Nm and Nn 2 are the attenuating hydrogen column 
densities, / is the partial covering factor, A is the normalization factor, -Efoid is the folding 
energy, V is the photon index, and cr(E) is the cross section of the neutral solar-abundance 
gas. It should be stressed that we adopt this model as a simple mathematical expression of 
the continuum spectra, and do not imply that it is the physical emission mechanism of the 
source. For example, it is possible that the continuum spectrum consists of more complicated 
components. The origin of the continuum emission will be discussed in forthcoming papers. 

We found that there is a complicated feature near 6 keV in the spectra of 1994 and 1995, 
which is not well accommodated by the continuum model alone. To fit the GIS continuum 
spectra in 1994 and 1995 satisfactorily, we had to exclude the energy range 6-8 keV. In 
the fit of the spectrum in 1996, the range 6-8 keV was included because no complicated 
feature was seen in there. The GIS spectra and the best-fit parameters are shown in Fig. 1 
and Table 2, respectively. As clearly seen, deep and broad absorption features exist at 6.5-7 
keV in the residuals of 1994 and 1995, as well as another absorption feature at ~ 8 keV. 
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In addition, there is a dip feature at 4 keV. In 1996, these absorption and dip features 
disappear or become shallower. 

To investigate these features, we use SIS spectra, which have a superior energy resolu- 
tion. The SIS spectra were fitted with the tentative continuum model determined from the 
GIS. Complicated absorption features were clearly seen in the spectra of 1994 and 1995 in 
Fig. 2. In particular, two narrow absorption lines (below and above 6.8 keV) are recognized 
in 1994. 



4.2. Absorption Line Features 

Next, we tried to fit the spectral features at the iron K energy band for the 1994 and 
1995 SIS data. The 5-10 keV range was initially fitted with a cut-off-power-law model, and 
Gaussians with negative normalizations were added one by one. Intrinsic Gaussian widths 
were allowed to be free but tied together. Four negative Gaussians were added to the models 
of 1994 and 1995, respectively, to achieve x 2 l v — 80.5/56 and 112.3/58, respectively, which 
are not satisfactory yet. We found that inclusion of an emission line at 6.4 keV to the 
model of 1994 and an absorption edge at 7.117 keV to that of 1995 improves the fits to an 
acceptable level. The best-fit parameters are shown in Table 3. The negative Gaussians are 
tentatively identified in Table 3 and referred to hereafter as negative Gaussians A, B, and 
so on. Adding an absorption edge to the model of 1994 or an emission line to that of 1995 
does not improve the fits. 

We also tried absorption edges instead of absorption lines, so that each negative Gaussian 
in Table 3 is replaced with an absorption edge. We found all the negative Gaussians cannot 
be replaced simultaneously, and at least one or two negative Gaussians are required to fit 
the spectra. These results suggest that absorption lines do in fact exist in the spectra of 
1994 and 1995. 

To check the consistency between the SIS and GIS, both data were fitted simultaneously. 
The equivalent widths and center energies of the emission and absorption lines, and optical 
depths of the edges were set equal for both data, while the parameters of the continuum 
model were allowed to be free. As a result, the best-fit parameters for the absorption features 
were consistent with those obtained from the SIS fit alone. Therefore, the absorption line 
parameters we have obtained are confirmed by both SIS and GIS, and we consider their 
detection to be unrelated to SIS pileup effects. 

To measure the absorption feature around the calcium K energy (~ 4 keV), the spectral 
continuum between 2.5 keV and 6 was fitted with an attenuated-power-law model. Then 
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an absorption-line model, negative Gaussian, was added to the fit model. It was found that 
the absorption feature of 1994 is fitted well with a negative Gaussian, while inclusion of a 
negative Gaussian does not improve the fit of 1995. The best-fit parameters are shown in 
Table 4. 



5. Discussion 

5.1. Interpretation of the Absorption Lines 

We notice the center energies of the negative Gaussians shown in Table 3 are consistent 
with or very close to the K-line energies of iron or nickel ions. Candidate resonant and 
inter-combination lines are listed in Table 5 (forbidden lines are not included). The energies 
of the negative Gaussians A and B coincide with those of Fe xxv Ka (ls 1 2p 1 1 P 1 and ls 1 2p 1 
3 Pi j2 in Table 5) and Fe xxvi Ka (2p x 2 Pi/2,3/2), respectively. The negative Gaussian C is 
considered to be Fe xxv K/3 (ls^p 1 ^i) or Ni xxvn Ka (ls l 2p l x Pi and ls^p 1 3 Pi, 2 ), 
although nickel is less abundant than iron in a plasma of cosmic abundance. The negative 
Gaussian D resides just between the Ni xxvm Ka (2p l 2 Pi/ 2 ,3/2) and Fe xxvi K{3 {Zp 1 
2 Pi/2) energies. These coincidences strongly support the interpretation that the absorption 
features consist of iron or nickel resonant absorption lines. Hereafter, we discuss the physics 
of the system assuming that the negative Gaussians A, B, C, and D are absorption lines of 
Fe xxv Ka, Fe xxvi Ka, Fe xxv Kf3 + Ni xxvn Ka, and Fe xxvi K/3 + Ni xxvm Ka, 
respectively, 

5.2. Curve of Growth Analysis 

The observed spectral features are obviously different from the P-Cygni profile which 
accompanies strong emission lines, and thus the absorption is probably not caused by a 
spherical wind. As for the origin of the absorption lines, we will consider line-scattering 
material anisotropically irradiated or anisotropically distributed around the central source. 
To estimate amount of the line-scattering matter in the line of sight, we calculate the "curve 
of growth", namely, the equivalent widths of an absorption line as function of the column 
density of scattering ions. 

We summarize here the basic equations we used to calculate the curves of growth. Let 
us assume particles with the column density N producing an absorption line at frequency 
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u by resonant scattering. The optical depth t{v) is expressed as 

t(u) = Ns<t> a (y\ (1) 
and the normalized equivalent width W v j uq is 

1 [ (1 - eM-T(v)])dv, (2) 
vo J 



Wu 



respectively, where s<p a {v) is the averaged cross section of scattering particles, normalized as 
J § a {y)dv = 1 (see Spitzer 1978, section 3.4). In CGS unit, the transition probability s is 
expressed as 

s = !L^! ( i_^exp[-^oAT]), (3) 
m e c bi 

where fi u is the oscillator strength, and the subscripts u and / denote values in the upper and 
lower states, respectively. The factor (b u /bi) exp[— hvo/kT], expressing the effect of induced 
emission, is defined to be equal to the actual rate of the density of ions in the upper state 
to that in the lower (Spitzer 1978, section 2.4). We assumed that the density of the ions in 
the upper state is negligible and omit the factor from the following calculation. The validity 
of the assumption is checked in § 5.6. If scattering particles obey the Maxwell-Boltzmann 
distribution with temperature T, (j>a{y) is written as a Voigt function H(a,u) with 

= vks H(a ' u) (4) 

H(a,«) = a - f Tl l 7 y2] X (5) 
7T J a 2 + (u- yy 

a - 4^I> (6) 

u = ^iy-vo) (7) 



b = J , 8 

V mi 

where A u i are the Einstein coefficients related to the oscillator strengths as 

(see Rybicki & Lightman 1979, section 10.6). Even if the energy distribution of particles is 
dominated by turbulence or other bulk motion rather than thermal motion, these equations 
are still applicable as a good approximation. In that case, T represents the velocity dispersion 
within the scattering gas. With these formulas and constants in Table 5, the equivalent 
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widths are calculated according to equation (2), and the curves of growth are plotted in 
Fig. 3. 

Using the curves of growth thus calculated, the column densities of the ions iV ion re- 
sponsible for the observed absorption lines may be obtained. Since Fe xxv Ka and Fe xxvi 
Ka lines are considered not contaminated by other ion species, iVpeXxv and iVpcXxvi are 
directly derived from the observed equivalent widths, assuming a temperature. From these 
column densities, equivalent widths of Fe xxv Kf3 and Fe xxvi Kf3 are estimated using the 
curves of growth. Then, these estimated equivalent widths are subtracted from the observed 
equivalent widths of the blends of Ni xxvn Ka + Fe xxv K/3, or of Ni xxvm Ka + Fe xxvi 
K[3, to obtain the equivalent widths of Ni xxvn Ka or Ni xxvm Ka. Column densities of 
N m xxvii an d N m xxvm are obtained using the curves of growth again. Ion column densities 
thus derived are shown in Table 6 for two extreme temperatures, 1000 keV and 0.1 keV. 

If the temperature is as high as 1000 keV (high-temperature limit), all the observed 
values of the equivalent width will be found in the "linear part" of the curves of growth, 
while if temperature is as low as 0.1 keV (low-temperature limit), they will be in the "square- 
root section" of the curves. It also should be noted that the column densities shown in the 
Table 6 are lower limits for both temperatures, since line photons scattered from matter 
located out of the line of sight would reduce the absorption line equivalent widths. That 
effect is not taken into account here. 

The derived column densities for the low-temperature limit are larger than those for the 
high-temperature limit by several orders of magnitude. If the column densities are so large 
as the low-temperature calculation suggests, and if abundances of the absorbing plasma is 
not very different from the solar value, the optical depth of the plasma for the Thomson 
scattering will exceed unity. This is unlikely, since the modification of the absorption line 
features that would be present due to Thomson scattering is not. It is plausible that either 
iron and nickel are significantly over- abundant, or that the kinematic temperature is as high 
as the high-temperature limit suggests. However, at such high temperature, iron and nickel 
would be fully ionized and absorption lines should not be observed. It may be that the 
ionization temperature of the absorbing plasma is lower than the kinematic temperature, 
or that the absorbing matter consists of several parts with different bulk velocities whose 
dispersion is comparable to thermal velocity of iron atoms at kT > 100 keV. In following 
discussion, we first consider the high-temperature limit in Table 6, which presumably gives 
the lower limit of column density. 
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5.3. Photo-ionized Plasma Model 

We note that the spectral parameters changed between the observations of 1994 and 
1995. The luminosity in the 1995 observation was higher than in 1994, and the ratio 
-^Fe xxvi /Np e xxv became larger, suggesting that the absorbing matter was more ionized. 
This behavior is consistent with characteristics of a photo-ionized plasma. Thus, we assume 
that the source of absorption lines in 1994 and 1995 is a photo- ionized plasma irradiated by 
the central source. In 1996, the system was considered to be in a different state, because, not 
only were absorption lines absent, but also the photon index and folding energy of the spec- 
tral continuum were different. We calculated electron temperature and ion-population dis- 
tribution of the photo-ionized plasma with XSTAR (Kallman & Krolik 1996), and searched 
for parameters which are consistent with those in the high-temperature limit in Table 6. 
To simplify the calculation, the plasma geometry was assumed to be isotropic around the 
source. It should be noted that assumption of isotropy is adopted only in the calculation, 
and that the realistic geometry of the line-absorbing plasma must be anisotropic, since oth- 
erwise absorption lines would be canceled by emission lines from out of the line of sight. In 
the calculations, abundance was assumed to be solar. The absorption-corrected 2-10 keV 
luminosity was estimated from Table 2 assuming a distance of 12.5 kpc, and used as total 
luminosity of the model. Continuum spectral parameters were drawn from Table 2. We 
arbitrarily fixed the ionization parameter £ = L/(nr 2 ) at the innermost boundary of the 
plasma to be 10 5 5 . Our results depend only weakly on this value, because iron and nickel 
is fully ionized with such a high ionization parameter at the innermost boundary, and thus 
column density of hydrogen-like or helium-like ions is not sensitive to the innermost condi- 
tion. The assumed geometry is schematically shown in Fig. 4. We searched a combination 
of plasma density and outer radius to give iV Fe xxvi an d N FeX xv consistent to observation. 
The determined parameters are shown in Table 7, and expected column density of ions are 
shown in Table 8. Obviously, the expected column density of nickel ions in Table 8 is much 
lower than the observation by one or two orders of magnitude. Therefore, we conclude that 
the relative abundance of nickel to iron is larger than that of solar value. 

We next searched for a combination of plasma parameters to reproduce the column 
density in the low-temperature limit in Table 6. Since the metal abundance was found 
to be non-solar, there was no reason to fix them to the solar value. We found that if 
both the iron and nickel abundances are multiplied by ~ 100, the resultant photo-ionized 
plasma gives the column densities in the low-temperature limit. Other parameters, such as 
the electron density, the ionization parameter at the outer boundary, and the outer radius 
of the plasma were hardly affected. Although the column densities are different between 
the high-temperature-limit and low-temperature-limit cases, the column-density ratios of Fe 
xxv to Fe xxvi are not much different between the two cases, and thus both yield similar 
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values of £1. (Subscript "0" and "1" denote values at the innermost and outermost boundary, 
respectively.) Thus we consider that £i and r\ in Table 7 are reliable, regardless of the 
kinematic temperature, the metal abundances, the electron density, or the innermost radius 
of the plasma. 

5.4. Geometry 

To produce absorption lines, the solid angle of the line-scattering material from the 
source must be small, otherwise the absorption lines would be offset by an emission line 
counterpart. From a moderate assumption that the solid angle is less than 47rexp[— 1], the 
disk half thickness, h, at the radius r\ is constrained as hjr\ < 0.4, where hjr\ is related to 
the solid angle as Q — 1/ \J\ + (h/ri) 2 . Therefore, the half opening angle of the plasma, 6, 
is constrained to be 9 < 20°. The inclination of the jets is 70° to our line of sight (Mirabel 
& Rodoriguez 1993). We might reasonably assume that the immediate neighborhood of the 
central source (e.g., the accretion disk) is also so inclined. From these two constraints, we 
conclude that the half thickness of the reprocessing plasma is ~ 0Ar±. Thus, total mass of 
the reprocessing plasma is estimated as ~ 4.4 x 10~ n M (1994) and ~ 1.5 x 10~ 10 M 
(1995). 

Note that n p in Table 7 is dependent on assumed geometry and has a large uncertainty, 
while £i and r\ are rather model independent and thus reliable. For example, a toroidal 
plasma distribution around the source can also explain the observed properties, and in this 
case the path in the plasma along the line of sight will be shorter. Therefore the proton 
density n p will be larger, conserving column density (ri — r ) x n p . The other parameters 
£i and r\ would not be changed much, because the former is determined from the ratio 
Np e xxvi /Np e xxv, and the latter is determined from the relation £i = L/{N c r\). 

5.5. Nature of the Reprocessing Plasma 

From the discussion above, we can construct a picture of the plasma producing the 
absorption lines. It is located at ~ 10 11 cm from the source, constrained to be within 9 < 20° 
of the disk plane, and has a total mass of 10~ n to 10~ 10 M . The kinetic temperature of 
iron ion is either high or low. If the metal abundance is comparable to the solar value, then 
the kinematic temperature must be as high as or higher than 100 keV, which corresponds 
to a velocity of ~ 10 s cm s" 1 . If the reprocessing plasma is hotter than 100 keV, it must be 
observed within 10 s of its production to avoid full ionization of iron and nickel. A post-shock 
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flow within 10 s after the shock is a candidate for such a hot plasma. Both the thermal motion 
of each ion and the bulk motion of the plasma may give such a high kinematic temperature. 
An inward or outward bulk flow with a velocity gradient or randomly moving blobs may be 
the X-ray reprocessor. 

As a candidate of such a flow, we suggest a radiation-driven disk wind. If gas at rest 
exists in the vicinity of the super-Eddington X-ray source, it would be blown out by Thomson 
scattering. The terminal velocity iWd would be 



where L is the luminosity of the source, orh is the cross section of Thomson scattering, r is 
the initial distance of the gas from the source, and m p is the mass of a proton. If such a gas is 
supplied from the accretion disk, we would observe a stream with a radial-velocity gradient 
from to fwmd- Substituting r = 10 11 cm and L = 4 x 10 38 erg s" 1 , a terminal velocity of 
^wind = 10 8 cm s" 1 is obtained, which coincides with the velocity inferred in the curve-of- 
growth discussion. It should be noted that most of the centroid energies of the absorption 
lines in Table 3 are blue-shifted from those in Table 5, which would be expected from the 
absorption lines from an outward flow with a velocity of 10 8 cm s _1 . Thus a radiation-driven 
disk wind at ~ 10 11 cm can well explain the observed absorption line features. Because the 
disk wind would be refreshed in every ~ 10 3 s, the outflow rate is estimated as 1CT 7 — 1CT 6 



On the other hand, if the plasma is very metal rich, the assumption of high-temperature 
is no longer necessary. The high column density of ions may be explained by low-temperature, 
metal rich plasma, if the metal abundance is 100 times higher than the solar value. Even in 
that case, the conclusion on the geometry of the plasma would not be changed appreciably. 
The nickel abundance of the jet material of SS 433 has been found to be a few tens times 
higher than the solar value (Kotani et al. 1997b), which suggests that anomalous abundances 
may not be rare in binary systems with strong jets. Thus they may be possible that iron 
and nickel are abundant in the environment of GRS 1915+105, and that the reprocessing 
plasma has a velocity dispersion smaller than 10 s cm s _1 or a lower temperature than 100 




(10) 



M yr 



-i 



keV. 



5.6. 



Alternative situations 



In the analysis above, the term of the induced emission (b u /bi) exp[— hvo/kT] in equa- 
tion (5.2) was neglected. The validity of the assumption can be confirmed as follows. Suppose 
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that a reprocessing plasma is illuminated by an external photon source. An ion in the plasma 
receives photons with a mean interval of l/(<j res F(i/ )), and emits each photon after a time 
~ l/A u i, where -F^o) is the photon flux of the source and <7 rcs is the resonant-scattering 
cross section. Therefore the ratio of the number density of ions at upper state to those at 
lower state is a res F{vo) / A u i. This is an overestimation because photons at vq are reflected 
by the plasma surface due to resonant scattering. If the ratio is much less than unity, the 
term of induced emission is negligible. Substituting a res = c 2 /(87rz/ 2 ), the ratio is estimated 
as 

which is much less than unity. Thus induced emission is shown to be negligible in the case 
of the photo-ionized plasma treated here. 

There is another case where the induced-emission term is important. The ratio ap- 
proaches unity on the surface of a neutron star (R = 10 6 cm). Under the presence of the 
induced emission, more ions are necessary to produce an absorption line. As an estimation of 
column density of Fe xxvi, we adopt 10 20 4 cm -2 from Table 6. The corresponding hydrogen 
column density would be 10 25 cm -2 , assuming solar abundance. To avoid full ionization of 
iron, the ionization parameter must be £ < 10 4 , and thus the density of the plasma must 
be > 10 22 cm" 3 , and the height of the atmosphere must be ~ 10 3 cm. Although the scale 
height of the atmosphere at kT ~ 100 is larger than 10 3 cm by an order of magnitude, such 
an atmosphere may produce absorption lines. However, the absorption lines made by the 
atmosphere would be gravitationally red-shifted, which was not observed. The atmosphere 
of a neutron star is unlikely to be the origin of the absorption lines. Thus we conclude that 
induced emission can be neglected in the usual situations where equation (5.6) is applicable. 



6. Summary 

Absorption lines of Ca XX Ka, Fe xxv Ka, Fe xxvi Ka, a blend of Ni xxvn Ka 
+ Fe xxv K/3, and a blend of Ni xxvm Ka + Fe xxvi K/3 have been discovered in the 
energy spectra of the transient jet source GRS 1915+105. These features can be explained 
by resonant scattering in a disk-like photo-ionized plasma. From the "curve-of-growth" 
analysis, column densities of the ions were determined as in Table 6. Adopting the values of 
high-temperature limit, the physical parameters of the plasma were determined to be: 

• outer radius T\ ~ 10 11 cm 

• half thickness h = 0.4ri 
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• density n p ~ 10 11 cm -3 

• ionization parameter at the outer boundary of the plasma £ = 10 3 ' 8 

• kinematic temperature of scattering ions > 100 keV, i.e., > 10 s cm s -1 , or metal 
abundance is as large as 100 times of the solar value. 

The kinematic temperature of the scattering ions should be very high to account for the 
observed large equivalent widths of the absorption lines, if solar abundance is assumed. 
If the temperature is so high, the plasma may not be in equilibrium, with an ionization 
temperature lower than the kinematic temperature. Thus the high kinematic temperature 
may be explained by a bulk motion, such that a velocity difference up to ~ 10 8 cm s -1 within 
the plasma flow, or perhaps plasma blobs moving randomly. We found that a radiation- 
driven accretion flow with a mass outflow rate of 10" 7 - 10" 6 M yr 1 can account for the 
absorption lines. Alternatively, if the plasma is extremely metal rich, such a high kinematic 
temperature is not necessary. Regardless of the kinematic temperature, we concluded from 
the observed iron to nickel abundance ratios that the abundances are non-solar. 



We thank Craig Markwardt for careful reading of the manuscript. 

Appendix A. Pile-up correction 

The SIS/ASCA suffers pileup when observing a bright source ~ 100 counts s _1 . If two 
or more photons fall on a single pixel or neighboring pixels (Moore neighbor) on a CCD chip 
of the SIS within the same readout period, they are recognized as a single photon by the 
on-board event-detection algorithm, and their energies are summed to make a single pulse 
height. This phenomenon is called "pileup" and causes a decrease in the photon counting rate 
and hardening of the spectrum. In observations of bright sources, pileup is non-negligible 
and correction is necessary for precise spectral analysis. We have developed a new code 1 for 
this purpose. This code also can cope with attitude fluctuations of the satellite, which were 
in fact significant during the GRS 1915+105 observation in 1994. 

Our code simulates the pileup using the observed events as the seeds. First, the source 
position is determined on the CCD image without correction, and the image core where 
the counting rate is higher than a certain threshold is discarded from further analysis. The 



Included in FTOOLS v5.0 available at http://heasarc.gsfc.nasa.gov/docs/software/ftools 
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threshold was empirically determined to be 0.023 counts pixel" 1 readout" 1 , above which we 
found the pileup effect was so significant that uncorrectable. The pixels on the CCD chip 
are divided into several annuli centered at the source. The annuli widths are determined 
so that the counting rate per pixel within each annulus is almost uniform. The number of 
detected events in the i-th annulus and the j-th readout is defined as a^. By dividing the 
chips into annuli, both the count rate and the image are eventually recovered. 

Next, the temporal order of all the iVj = events in the i-th annulus is randomized, 
while each event retains its position, pulse-height, and charge-distribution pattern, i.e., grade 
information. The event series before and after randomization are denoted as 

7a, 7i2, • • • ,7iNi (12) 

and 

7w7(l), 7«t(2), ••• ,7i<r(Ni), (13) 

where cr(k) is a randomizing function which returns an integer between 1 and iVj, and satisfies 
u{k) a (I) if k ^ I. For each read-out, the randomized events are injected on a virtual 
CCD chip, and the same event detection algorithm as that on-board is applied (Gendreau 
1995). For example, to simulate the j-th readout of the i-th ring, a series of events 

is created (note that Y? k=1 dik — S^Ojfe = a^), and charge-distribution patterns on the virtual 
CCD chip are calculated from the position, pulse-height, and original charge-distribution 
pattern or grade of each input event. If the SIS observation mode is BRIGHT and the 
information on split charge has been lost, the charge is assumed to be 40 ADU (the split- 
threshold). The spectrum of the injected events is recorded as /2(E). If there is no pileup, 
all dij events should be detected. On the other hand, the detected events may be less than 
Ojj if pileup is accounted for. In that case additional events are randomly picked up from the 
total Ni events and injected onto the virtual chip (also added to /2(E)) until events are 
detected. When the number of detected events has reached a^, the spectrum determined 
by the event detection algorithm is recorded as /3(E), which will be the simulated spectrum 
affected by the pileup. This procedure is repeated for all the annuli and readouts. It should 
be noted that all the original events are used in the whole process at least once. By using all 
the events, fluctuations caused by random pick up will be reduced to minimum. Finally, the 
pure pileup spectrum fi{E) and resultant pileup-corrected spectrum f*,{E) are estimated as 

h = /s-/ 2 (15) 

h = /i-/4 = /i + / 2 -/ 3 , (16) 

where fi(E) is the observed spectrum before the correction. In Fig. 5, an example of the 
pileup correction is exhibited. 
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Table 1: Observation log 



No. 


Start 


End 


Exposure Time a 


SIS Mode 




(UT) 


(UT) 


(ks) 




1 


1994/09/27 00:40 


1994/09/27 13:00 


10 


BRIGHT 


2 


1995/04/20 09:50 


1995/04/20 22:40 


16 


FAINT 


3 


1996/10/23 02:40 


1996/10/23 18:20 


11 


BRIGHT 


4 


1997/04/25 08:10 


1997/04/25 21:00 


20 


BRIGHT 


5 


1998/04/04 16:20 


1998/04/05 03:30 


20 


BRIGHT 


6 


1999/04/15 20:10 


1999/04/16 15:00 


20 


BRIGHT 



a Bit HIGH only. 



Table 2: Best-fit parameters of the GIS spectral continuum 



No. 


N m 
(10 22 cm- 2 ) 


/ 


^H2 

(10 22 cm" 2 ) 




-Efold 

(keV) 


r 


xl iy) 


l b 


2 Q7+ - 10 
z - y ' -0.10 


n KK 1+0.040 
U.OOl_ 040 




4 fi7+ - 44 


1 no+0.04 
i - zo -0.04 


-0.911H 


1.13 (131) 


2 b 
3 


3 46 +0 ' 10 
°-^ u -o.io 

4 - ' °-0.04 


U.OOO_ 070 

n 2 -| q+0.066 
u.z,±y_ 045 


81^ 

71 +22 

1 r_i4 


5.8418;" 


i qc+0.05 
i -° d -0.05 

4 61 +0 - 28 


-1-18^ 

i - ou -0.08 


1.13 (131) 
1.12 (174) 



a In unit of photons s 1 cm 2 keV 1 at 1 keV. 
h The Energy range 6-8 keV was excluded. 
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Table 3: Best-fit parameters of the feature around iron-K absorption 



Model Component 


Parameter 


1994 


1995 


Emission Line 


Energy (keV) 


6.4841°;°^ 






Width (eV) 


(fixed) 






EW (eV) 


15.8123 




Negative Gaussian A 


Energy (keV) 


6.7201°;^ 


6.6901°;°^ 




Width a (eV) 


< 39 


< 89 




EW (eV) 




27-5lg;3 


Negative Gaussian B 


Energy (keV) 


7.0001°;°?° 


6.9641°;°- 




EW (eV) 


43.2l|| 




Negative Gaussian C 


Energy (keV) 


7 843+ - 037 


7 oon+0.043 

i -ooy_ 037 




EW (eV) 


30.8+^g 2 


29.9ll;| 


Negative Gaussian D 


Energy (keV) 


8-i99±8:8S 


8 181 +0.060 




EW (eV) 




1Q+9 
iy -10 


Absorption Edge 


Energy (keV) 




7.117 (fixed) 




Max. Opt. Depth 




n i 99+O.O27 
U.±ZZ__q Q27 


xl M 




1.03 (56) 


1.05 (57) 



"Upper limit at 1 a confidence level. Set free but same for all negative Gaussians. 



Table 4: Best-fit parameters of the absorption features at 4 keV 



Model Component 


Parameter 


1994 


1995 


Negative Gaussian X 


Energy (keV) 


4 n7 i +0.068 
4 - u ' i -0.097 


4.05l°- 




Width 


(eV) 


< 210 


10 (fixed) 




EW 


(eV) 


10.2115 


< 7.9 


xl H 






0.80 (56) 


1.24 (57) 
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Table 5: Absorption-line candidates 



Ion 


Upper Level 


Energy 




Oscillator Strength 


Einstein Coefficient 






(keV) 






(s" 1 ) 


Fe xxv 


ls^p 1 3 Pi 


6.6684925 (1) 




0.0687 (1) 


4.42 x 10 13 


Fe xxv 


ls^p 1 1 P 1 


6.7011266 (1) 




0.703 (1) 


4.57 x 10 14 


Fe xxvi 


2P 1 2 Pl/2,3/2 


6.9519639, 6.9731781 


(3) 


0.44 (4) 


1.5 x 10 14 


Ni xxvii 


ls l 2p l 3 Px 


7.7668938 (1) 




0.0883 (1) 


7.70 x 10 13 


Ni xxvii 


ls^p 1 1 P 1 


7.8062340 (1) 




0.683 (1) 


6.02 x 10 14 


Fe xxv 


ls^p 1 x Pi 


7.8820244 (2) 




0.138 (2) 


1.24 x 10 14 


Ni xxviii 


2P 1 2 Pl/2,3/2 


8.0731039, 8.1017429 


(3) 


0.44 (4) 


2.07 x 10 14 


Fe xxvi 


3P 1 2 Pl/2,3/2 


8.2636944, 8.2526875 


(3) 


0.046 (4) 


2.3 x 10 13 



References. — (1) Drake 1979; (2) Lin, Johnson, & Dalgarno 1977; (3) Sugar & Corliss 1985; (4) Mewe, 
Groncnschild, & van den Oord 1985, and references therein. 



Table 6: Column density of absorbing ions 





High- Temperature Limit (kT - 


= 1000 keV) 


1994 
1995 


Fe xxv Fe xxvi Ni xxvii 

17.7+^ 18.018:1 i7.5i8;| 
17.618:1 i8.H8:l 17.5+8J 


Ni xxviii 

17.6I8:? 

1 7 4+0-5 




Low- Temperature Limit (kT 


= 0.1 keV) 


1994 
1995 


Fe xxv Fe xxvi Ni xxvii 
19.3+81 20.3+81 18.8+81 

19.0I8I 20.4I8I i8.8l8i 


Ni xxviii 

19.4+°J 

i9.il?:! 



Note. — In units of logio[cm 2 ]. 
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Table 7: Parameters of the photo-ionized plasma 





L 






r a 


6 


ri a 


n p 




(erg s 




(erg s _1 cm) 


(cm) 


(erg s _1 cm) 


(cm) 


(cm -3 ) 


1994 


4.0 x 10 38 


(fixed) 


10 5 - 5 (fixed) 


9 x 10 10 


10 3.8 


6 x 10 11 


1.5 x 10 11 


1995 


8.8 x 10 38 


(fixed) 


10 5 - 5 (fixed) 


14 x 10 10 


1Q 3.8 


9 x 10 11 


1.5 x 10 11 



Calculated from ionization parameter £, and not a free parameter. 



Note. — Subscript "0" and "1" denote values at the innermost and outermost boundary, respectively. 



Table 8: Model prediction of ion column density 





H II 


vin 


Si xiv S xvi 


Ca xx 


Fe xxv 


Fe xxvi 


Ni xxvii 


Ni xxviii 


1994 


22.9 


16.1 


15.9 16.0 


15.9 


17.7 


18.0 


15.1 


16.7 


1995 


23.1 


16.4 


16.0 16.1 


15.9 


17.6 


18.1 


16.7 


16.9 



Note. — In unit of logio[cm~ 2 ]. Ion species less than 10 16 cm" 2 were omitted, as were naked ions. The 
abundance was assumed to be solar. Input parameters were determined so that Fe xxv and Fe xxvi were 
consistent with the observations. 
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Energy (keV) 



(c) 1996 October 



Fig. I.— GIS Spectra, (a) 1994 September, (b) 1995 April, and (c) 1996 October. Observed 
data (cross), best-fit model (histogram), and residuals (bottom panel) are plotted. Each 
error bar includes systematic uncertainties of 2 %. The data of GIS-2 and GIS-3 were 
combined. Dead time and telemetry saturation were not corrected. The energy range 6-8 
keV in the data of 1994 and 1995 was excluded in fitting. The structure at 4.7 keV is due 
to an instrumental Xe-L edge. 
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(a) 1994 September 



(b) 1995 April 



(c) 1996 October 




1 1 Z 5 10 

Energy (keV) 



Fig. 2.— SIS Spectra, (a) 1994 September, (b) 1995 April, and (c) 1996 October. Observed 
data (cross), best-fit model (histogram), residuals (crosses in lower panel), and residuals 
for the model without absorption lines (solid line in lower panel) are plotted. Absorption 
features are clearly seen in the solid curve. As for the data of 1994 and 1995, only the 
range around the iron absorption features are plotted, while the overall spectrum is plotted 
for 1996. The data of SIS-0 and SIS-1 were combined. The structure around 2-2.5 keV is 
instrumental. 




16 18 20 22 16 IB 20 22 

lo g N Fe XXV ( cm_2 ) N Fe XXVI ( cm ~ 2 ) 

(a) Fe xxv Ka (b) Fe xxvi Ka 




16 18 20 22 16 18 20 

lo g N M XXVII < Cm_Z ) lo S N M XXVIII ( cm_2 ) 

(e) Ni xxvii Ka (f) Ni xxvm Ka 



Fig. 3. — Curves of growth, (a) Fe xxv Ka, (b) Fe xxvi Ka, (c) Fe xxv Kf3, (d) Fe 
xxvi K(3, (e) Ni xxvii Ka, and (f) Ni xxvm Ka. The curves in each panel correspond to 
kinematic temperature of 1000 keV, 100 keV, 10 keV, 1 keV, and 0.1 keV, from high to low. 
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Fig. 4. — Assumed geometry. An X-ray source (gray sphere) is embodied in a plasma disk 
expressed as a frame. 
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Fig. 5. — Spectrum before and after pile-up correction. The Sensor-0 Chip-1 data of 1994 
are presented as an example. The observed spectrum before correction fi(E), the calculated 
pure pileup component f±{E), and the resultant pileup-corrected spectrum f${E) are shown. 
Energy is expressed in analog-to-digital-converter units (ADU), and the range is between 
0.5-10 keV. In the lower energy band (50-200 ADU), photons are lost by pileup and thus 
fi{E) is negative by definition. In the higher energy band (> 250 ADU), pseudo photons 
are made by pileup. Although it is difficult to discern by eye in a linear plot, the effect in 
the higher band is not negligible. 



-25 - 



Erratum: "ASCA Observations of the Absorption-Line Features 
from the Super-Luminal Jet Source GRS 1915+105" (ApJ, 539, 

413 [2000]) 

We have discovered errors in the calculation of some of the Einstein coefficients in Table 5 
and some plots in Figure 3. Due to the errors, the square-root section of the curves of growth 
of Fe xxvi Ka and Kf3, and Ni xxvm Ka in Figure 3 were a few times underestimated. We 
correct Table 5 and re-plot Figure 3 together with the unaffected curves 2 . We thank Aya 
Kubota and Chirs Done for pointing out this error. 



2 The program to calculate the corrected curves of growth is available at 
http : / /www . hp . phys . titech . ac . jp/kotani/cog/index . html 
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TABLE 5 
Absorption-line candidates 



Ion 


Upper Level 


Energy 


Oscillator Strength 


Einstein Coefficient 






(keV) 




(s" 1 ) 


Fe xxv 


ls l 2p l 3 Pi 


6.6684925 (1) 


0.0687 (1) 


4.42 x 10 13 


Fe xxv 


ls^p 1 


6.7011266 (1) 


0.703 (1) 


4.57 x 10 14 


Fe xxvi 


2P 1 2 Pi/2, 3 /2 


6.9519639, 6.9731781 (3) 


0.44 (4) 


3.1 x 10 14a 


Ni xxvii 


ls l 2p l 3 Px 


7.7668938 (1) 


0.0883 (1) 


7.70 x 10 13 


Ni xxvii 


ls^p 1 1 P 1 


7.8062340 (1) 


0.683 (1) 


6.02 x 10 14 


Fe xxv 


Is^p 1 x Pi 


7.8820244 (2) 


0.138 (2) 


1.24 x 10 14 


Ni xxvin 


2P 1 2 Pl/2,3/2 


8.0731039, 8.1017429 (3) 


0.44 (4) 


4.1 x 10 14a 


Fe xxvi 


3P 1 2 Pl/2,3/2 


8.24636944\ 8.2526875 (3) 


0.046 (4) 


4.5 x 10 13a 



"Corrected. 
b Typo corrected. 

References. — (1) Drake, G. W. F. Phys. Rev. A, 19, 1387 [1979]; (2) Lin, C. D., Johnson, W. R., & 
Dalgarno, A. ApJ, 217, 1011 [1977]; (3) Sugar, J., & Corliss, C. J. Phys. Chem. Ref. Data, 14, Supplement 
No. 2 [1985]; (4) Mewe, R., Gronenschild, E. H. B. M., & van den Oord, G. H. J. A&AS, 62, 197 [1985], and 
references therein 
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16 18 20 22 16 18 20 22 

l °S N Nl xxvn (cm- 2 ) log N m xxvm (cm- 2 ) 

Fig. 3e Fig. 3f 



Fig. 3 — Curves of growth, (a) Fe xxv Ka, (b) Fe xxvi Ka, (c) Fe xxv K/3, (d) Fe xxvi K(3, 
(e) Ni xxvn Ka, and (f) Ni xxvm Ka. The curves in each panel correspond to kinematic 
temperature of 1000 keV, 100 keV, 10 keV, 1 keV, and 0.1 keV, from high to low. 
This figure replaces Fig. 3 of the original manuscript. The square- root section of the curves 
in Fig. 3b, 3d, and 3f are a few times increased by this correction. The unaffected figures, 
3a, 3c, and 3e are also shown. 



